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Abstract 

The distribution of the QCD topological charge can be described by cumulants, with the 
lowest one being the topological susceptibility. The vacuum energy density in a 0-vacuum 
is the generating function for these cumulants. In this paper, we derive the vacuum energy 
density in SU(2) chiral perturbation theory up to next-to-leading order keeping different up 
and down quark masses, which can be used to calculate any cumulant of the topological charge 
distribution. We also give the expression for the case of SU(A) with degenerate quark masses. 
In this case, all cumulants depend on the same linear combination of low-energy constants 
and chiral logarithm, and thus there are sum rules between the A-flavor quark condensate 
and the cumulants free of next-to-leading order corrections. 
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1 Introduction 


Because of the axial U(l) anomaly, there exists a 0-term in quantum chromodynamics (QCD) 
which is a topological term. The partition function of QCD in a 0-vacuum is given by 

Z{e)= ( 1 ) 

where ^gcoiG, q, g] is the QCD action at 0 = 0 with G and q being the gluon and quark fields, 
respectively, and Q is the topological charge 

Q = / d^xG^’'ix)GP^x) , (2) 

with G^''{x) the gluon field strength tensor. In the Euclidean space with a finite space-time volume 
E, the partition function Z{6) is dominated by the ground state, i.e. vacuum, energy of QCD for 
large enough V (see, e.g. Ref. [1]), and we have 

Z{0)=e-^^-^^\ or evac(0) = -^ln^(0), (3) 

where CvaciO) is the vacuum energy density in the 0-vacuum. The distribution of the topological 
charge can be described in terms of moments, which are the expectation values (Q^")g_Q with 
positive integer n, or cumulants defined as 


C2n — 


d^”evac(fi') 


9=0 


(4) 


The leading cumulant is the topological susceptibility, C 2 = Xt- It and the fourth cumulant are 
given by the well-known formulae 

X- = VW'>.-. = (5) 


These topological quantities are important to understand the QCD vacuum as well as to extract 
physical observables from lattice simulations at a hxed topology [1, 2]. They can be measured on 
lattice using various methods, see, e.g., Refs. [3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19]. 

For large volume and small quark masses, the strong interaction dynamics is determined by the 
Goldstone bosons originating from the spontaneous breaking of the light-quark chiral symmetry, 
and thus can be well described by chiral perturbation theory (CHPT) [20, 21]. Both of Xt and C4 
have been calculated in CHPT in both leading order (LO) and next-to-leading order (NLO) [22, 
23, 1, 24, 25, 26, 27, 28]. Earlier discussions in the large Nc limit can be found in Refs. [29, 30]. 
The NLO calculations for xt in Refs. [25, 28] and for C4 in Ref. [28] were performed for an arbitrary 
number of flavors with different masses, and based on the generating functionals of CHPT [21] 
expanded around 0 = 0 up to 2-point loops (up to 1-point tadpole loops for the topological 
susceptibility [25]). 

In this paper, we will derive a general formula for the vacuum energy density in SU(2) chiral 
perturbation theory keeping different masses for the up and down quarks. The derivation involves a 
direct calculation of the logarithm of the determinant for the free Goldstone bosons in a 0-vacuum, 
and thus does not require an expansion up to a finite n-point loops. In this sense, it contains a 
summation of all one-loop diagrams at NLO in the chiral expansion, i.e. O (p^) with p denoting 
a small momentum or Goldstone boson mass, contributing to the vacuum energy. The expression 
for the vacuum energy density can then be used to calculate any cumulant of the distribution of 
the QGD topological charge defined in Eq. (4). 

It was emphasized in Ref. [28] that lattice simulations of these topological quantities with 
degenerate quarks are very interesting to pin down the N"-flavor quark condensate. Although both 
the topological susceptibility and the fourth cumulant depend on several low-energy constants 
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(LECs) in the NLO chiral Lagrangian, in addition to the quark condensate, the authors found 
an interesting linear combination, xt + iV^C4/4 with N the number of flavors, independent of any 
LEG. Thus, such a combination is particularly suitable for extracting the iV-flavor averaged quark 
condensate whose absolute value is 

Sat = F^Bn , (6) 

where F^, the pion decay constant, and Bat are defined in the chiral limit. For determinations 
of the quark condensate from lattice calculations of the topological susceptibility, we refer to 
Ref. [26, 16, 17]. Stimulated by this insight, we will also derive general expressions for the SU(IV) 
vacuum energy density and cumulants with degenerate quarks. It turns out that all the cumulants 
depend on the same linear combination of the NLO LECs and chiral logarithm. As a result, one 
can construct linear combinations of the cumulants free of NLO corrections. 

At this point, we notice that higher cumulants can be obtained from lower ones and moments 
using the following recursion relation 


C2n — ( 1) 


n+1 




n —1 


V 




771—1 


2n-l 
2m — 1 




n—m) 


C2m 


J 0=0 


(7) 


2 Vacuum energy in SU(2) chiral perturbation theory 

2.1 Leading order 

Because the d-angle can be rotated to the phase of the quark mass matrix by an axial U(l) rotation, 
the 0-dependence of physical quantities can be studied by using a complex quark mass matrix. At 
LO, O (p^), of SU(iV) chiral perturbation theory, the vacuum energy density in a 0-vacuum for N 
quarks is given by 

e^^M = -^{x0UUxlUo). ( 8 ) 

where xe = 2 BmM. exp(i0/A^) with A4 being the real and diagonal quark mass matrix, and the 
vacuum alignment Uq can be parametrized as a diagonal matrix Ug = diag{e®‘^L 
with the constraint (pi = 0. The angles (ft are determined by minimizing the vacuum energy. 
It is equivalent to removing the tree-level tadpole terms of the neutral Goldstone bosons which 
would induce vacuum instability [31, 32, 33]. 

In this section, we will study the case with N = 2. We will drop the subscripts in F 2 and B 2 
to be consistent with the traditional notation in GHPT. With Ug = diagle®"^, 6“*''’}, we have 


/o', 2 - f ^ ^ \ 

6vac(^) = 2E Bm ( cos - cos p — e sin - sin 1 


(9) 


where m = (m„-|-md)/2isthe average mass of the up and down quarks and e = (md — mu) / (mu + 
md) quantifies the strong isospin breaking. Minimizing the vacuum energy with respect to p, one 
gets [1] 


tan ip = —e tan — . 


( 10 ) 


Substituting this into Eq. (9), we get the vacuum energy density at LO, up to an additive nor¬ 
malization constant [1] 

e(i)(0) = -F2M2(0), (11) 


where M^(0) is the LO pion mass squared in a 0-vacuum 


[ 1 ] 


(0) = 2Bm cos - \ 1 + tan^ 


( 12 ) 
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Notice that in the absence of the electromagnetic interaction, the neutral and charged pions have 
the same mass at LO. The cumulants of the distribution of the topological charge can then be 
easily obtained. For instance, the topological susceptibility and the fourth cumulant at LO are 

= lF^Brh{l-e^), =-^-F^Bfh {l + , (13) 

which have been derived before in Refs. [22, 25, 24]. 


2.2 Next-to-leading order 

At NLO, there are contributions from both the tree-level terms in the O {p'^) chiral Lagrangian 
and one-loop diagrams. The vacuum energy density up to NLO is given by 

evac(0) = e(^) (0) + + e(if-)(0), (14) 

where e™l(0) is given in Eq. (9), (0) is the one-loop contribution to be calculated later on, 

and the NLO tree-level contribution is 


=("’t’'ee)(0) = -^(^xlUo+xeU^oY + ^ (xj Uo - X0 U, 


hi 


‘3 / t 
,X9 Xe 


hi — /i3 


Re(det xe) 


= -M\0){k + l7 


(1 — e^) tan(0/2) 


n 2 ' 


1 -I- tan^(0/2) 

-2R^to^ [(hi -I- hs) (l -I- e^) -I- (hi - ha) (l - e^) cos0] 


(15) 


where I 3 , Ir and hi, ha are the LECs and high-energy constants (HECs), respectively, in the NLO 
two-flavor chiral Lagrangian [20],^ and we have used Eq. (10).^ Because both I 3 and hi are 
ultraviolet (UV) divergent [20], 




hi = h\ + 2 \, 


(16) 


with A the divergence at the space-time dimension d = 4 in dimensional regularization, 

,d-4 


A = 




1071 ^ 




(17) 


where /r is the scale in dimensional regularization, eitc(0) is UV divergent as well, and the 
divergence is (the divergence is nontrivial in a 0-vacuum noticing the 0-dependence) 

3A 

g(4.tree.oo)(0)^__^4(^)^ (18) 

As will be shown, this divergence is exactly cancelled by the one from loops in elic°°^^(0). 

Before proceeding to calculating the loop contribution to the vacuum energy density, let us 
discuss the main difference between our treatment (see below) and the one in Refs. [25, 28]. In 
those papers the authors took the expression of the generational functional in Ref. [21]. It is 
normalized to the free fields at 0 = 0 (notice that Refs. [20, 21] assume 0 = 0). Then the loops 

^Here we use the SU(2) X SU(2) notation rather than the 0(4) one in the original paper, see, e.g., [34]. 

^In principle, the vacuum alignment determined by minimizing the LO vacuum energy gets shifted due to the 
presence of the higher order terms, Zy in this case. However, this shift only provides a perturbation and is of one 
order higher compared to the angle in Eq. (10). It introduces CP-odd vertices (see, e.g., Refs. [35, 36]) and 
does not affect CP-even quantities up to O (p^), thus irrelevant for us. It is for this reason that the topological 
susceptibility up to NLO in the chiral expansion calculated in Ref. [28] agrees with that in Ref. [25], where the 
vacuum alignment was calculated by minimizing the LO and NLO vacuum energy, respectively. 
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were calculated using the Goldstone boson masses at 0 = 0, and the 0-dependence is kept in the 
operator cr^ defined as (we have replaced U containing quantum fiuctuations of Goldstone bosons 
by C/q relevant for the vacuum energy) 

= 5 ({Ap, (xJ C/o + X« C/t) ) - SpqMUq) , (19) 

where \p are linear combinations of the S\J{N) generators introduced to diagonalize the LO mass 
term [21], and Mp(0) are the LO Goldstone boson masses at 0 = 0. This amounts to an expansion 
around 0 = 0, which is perfectly fine for the calculation of the cumulants of the topological charge. 
However, we notice that the first term in the above equation is in fact SpqMp^O). If we expand 
the one-loop generating functional around the one for the free fields in a 0-vacuum, 


ZoiO) = ^ Indetllo)^') = ^Trln£)o(0), 


( 20 ) 


where Tr stands for taking trace in both the flavor (this is the space of the adjoint representation 
which is 3-dimensional for the SU(2) case) and coordinate spaces, and Dq{6 ) is a differential 
operator. 


Dopq{ 0) = SpQ 


a^5^ + M2(0) , 


( 21 ) 


then Mp{0) in Eq. (19) needs to be replaced by Mp{9) and (Jpg vanishes. As a result, the only 
term left in the one-loop generating functional relevant for the vacuum energy is Zq{6 ). Thus, the 
vacuum energy density is given by 


eii>°P)(0) = -^Trlni?o(0). (22) 

For the case of SU(2), because the neutral and charged pions have the same mass at LO, 
Mp(0) is given by Eq. (12), and Do{0) = Isxs [d^d'^ + M^(0)], where the unit matrix has the 
dimension of the adjoint representation for SU(2). Extending these considerations to the case of 
N degenerate quark flavors and using dimensional regularization, we obtain 


(4,loop) 
’^vac 






-p^ + M^{ 0 ) 


(iV^ - 1)M^{9) 


X 

2 


1 

1287r2 


1 


— 2 In 



(23) 


For N = 2, one sees that the UV divergence cancels exactly the one in Eq. (18). The sum of 
Eqs. (9), (15) and (23) provides the vacuum energy density in a 0-vacuum up to NLO, 


evac(0) = -F^M\9)-M\9) 


1287r2 


1 - 2 In 


M2(0) 


+ /" + hC - ha + /t 


(1 — e^) tan(0/2) 


n 2 


1 -I- e2 tan^(0/2) 


(24) 

where we have dropped 0-independent constant terms. The renormalized LEG and HEC h\ are 
scale dependent [20] and this scale dependence cancels that in the chiral logarithm resulting in a 
scale-independent vacuum energy density in a 0-vacuum. This is the main result of our paper. It 
is then trivial to obtain the expression for any cumulant, and the lowest two are 


+ 2.^ - 34) + (1 - 4) (I - 4) + yb ^ 

-2 [/^ + hl-h3-h{l + 2e^ - 3e^)] 1 + 0 (p®) . (25) 
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They agree with the general 7V-flavor expressions in Ref. [28] for = 2. Furthermore, in the isospin 
symmetric case, they depend on the same combination of the LECs and HECs, — It + h\ — h^. 


3 SU(A/^) with degenerate quark masses 


The evaluation of the functional determinant Zq{6 ) or in Eq. (23) only requires the 

Goldstone bosons to be degenerate. Therefore, it is easy to generalize the result in the previous 
section to the case of SU(A^) with degenerate quark masses.^ The one-loop contribution to the 
vacuum energy density in a 0-vacuum is given by Eq. (23) as well with M{6) replaced by the LO 
Goldstone boson mass for SU(A^), see below. 

When all the quarks are degenerate with a mass m, the vacuum is given by Uq = Inxn- 
With the O (p^) Gasser-Leutwyler Lagrangian for SU(A^) [21], we get the tree-level contribution, 
including both LO and NLO, to the vacuum energy density in a 0-vacuum 


e 


N 


e 


26 


= —NF^Bnui cos — — A.NBj^rrr ( diVLe cos'^ — — ANL^ sin-h 2Lg cos --— I- 4iL2 , 


N 


N 


N 


(26) 

where Ley.8 are LECs and H 2 is a HEC. Among them, Lq,Lt and H 2 contain a UV divergent 
piece which can be calculated using the heat kernel method with path integral [21, 38] 




N^ + 2 


A, 




7V2 - 4 


A, 


H 2 = m 


-4 


A. 


(27) 


I 6 A 72 ° 

It is straightforward to check that these divergences cancel the one in eiac°°^^ in Eq. (23). The 
vacuum energy density in a 0-vacuum up to NLO is then 

evac(0) = 


f - 1 


/ 9 ^ \ 1 

1 1287r2 

FAN iNLlFLl-NLrtan^—] 



(28) 


with the scale-dependent finite LECs Lg and Lg, where M^{6) = 2BNmcos{9/N), and the cu- 
mulants are 


C2ri = 


(- 1 ) 


n+1 


N^n-l 

- 1 


F^Bjsim F 

Blf-m'^^N,2n 


- 1 
647r2lV 


1 - 2 In 


2 BMm\ 

) 


-f 8(A^ Lg -|- Lg -f N Lt) 


IGtt^ 

with the number ^N, 2 n defined as 

^N, 2 n = 


d02'’ 


COS 


2 0 1 / 0 
— In cos — 


N 


N 


(29) 


(30) 


0=0 


One sees that all cumulants depend on the same linear combination of the LECs, as observed in 
Ref. [28] for the topological susceptibility and the fourth cumulant, and chiral logarithms. From 
this it is easy to construct LEC-free combination of cumulants which can be used for a clean 
extraction of the N-flavor quark condensate from lattice simulations as suggested in Ref. [28]. 
Examples are 


Xt 


Xt 


7V2 

7V4 

l6 


C4 


■C6 


SF^NBNm 3 {N^ - 1) B%m^ 

AN 327r2iV2 +[P ) ^ 

ISE^Rjym 15 (A^2 _ 1) g 

16fV 647r2fV2 + IP j 


(31) 


^For SU(-/V) with different quark masses, one may expand around 0 = 0, \nDo(0) = InZ)o(O) + D^^(O) A(0) + 
Dq^{0) A(0) D^^(O) A(0) + ... with Apq{ 0) = SpQ — Mp(0)j. This gives the general formulation used in 

Ref. [28]. 
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where the first expression was already proposed in Ref. [28]"'’. More interestingly, the NLO correc¬ 
tions can be canceled out completely in certain linear combinations, and lead to sum rules between 
the QCD topological sector and the spontaneous breaking of chiral symmetry, such as 


Eat = 


N /8 


2N'^ 


7Xt 


-C4 


IE 


C6 


o(/). 


(32) 


In fact, in the chiral limit, we have the following exact relation as can be seen from Eq. (29) 

Sat = 7 rp(0) = lim (-1 )^+^, (33) 

m-^0 m 

where we have displayed the Banks-Casher relation [39] linking the quark condensate to the zero¬ 
mode spectral density of the Euclidean Dirac operator, denoted by p(0), as well. These relations 
can be simply obtained using the LO expression for the vacuum energy density, and suggest that 
there is an intimate link between the QCD topological sector and the spontaneous breaking of 
chiral symmetry. 


4 Summary 

We have derived the expressions for the vacuum energy density in a 0-vacuum in SU(2) CHPT up to 
NLO keeping different up and down quark masses as well as in SU(iV) CHPT with degenerate quark 
masses. They can be used to calculate the cumulants of the QCD topological charge distribution 
which are important quantities to study QCD in the low-energy strong coupling regime. In the 
case of degenerate quark masses, all cumulants depend on the same linear combination of low- 
energy constants, as already observed for the topological susceptibility and the fourth cumulant 
in Ref. [28]. Therefore, one can construct many combinations of the cumulants depending only 
on the quark mass and condensate. They can be used to extract the quark condensate in lattice 
simulations without contamination from LECs. Furthermore, we find sum rules relating the quark 
condensate to the cumulants free of NLO corrections. It would be interesting to check such 
relations in lattice QCD. 
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